arXiv:1507.01964vl [cond-mat.soft] 7 Jul 2015 


Fast and rewritable colloidal assembly via field synchronized particle swapping 

Pietro Tierno/’Tom H. Johansen,and Thomas M. Fischer^ 

Estructura i Constituents de la Materia, Universitat de Barcelona, Spain 
Institut de Nanociencia i Nanotecnologia IN^ UB, Universitat de Barcelona 
^''Department of Physics, University of Oslo,P. O. Box lOfS, Blindern, Norway 
Institute for Superconducting and Electronic Materials, University of Wollongong, Wollongong, NSW 2522, 
Australia 

Institut fiir Experimentalphysik V, Universitat Bayreuth, 95440 Bayreuth, Germany 
(Dated: 9 July 2015) 

We report a technique to realize reconfigurable colloidal crystals by using the controlled motion of particle 
defects above an externally modulated magnetic substrate. The transport particles is induced by applying a 
uniform rotating magnetic field to a ferrite garnet film characterized by a periodic lattice of magnetic bubbles. 
For filling factor larger than one colloid per bubble domain, the particle current arises from propagating defects 
where particles synchronously exchange their position when passing from one occupied domain to the next. 
The amplitude of an applied alternating magnetic field can be used to displace the excess particles via a 
swapping mechanism, or to mobilize the entire colloidal system at a predefined speed. 


The flow of electrical current in a conductor arises 
when electrons at the Fermi level are scattered from oc¬ 
cupied states into unoccupied ones due to the interaction 
with an external electric field. Following this electronic 
analogy, it is a compelling idea to find ways to separate 
an ensemble of identical colloidal particles into a set of 
immobile low energy particles, and colloids which can 
become mobile due to an external field. When applied 
to microtechnological devices like lab-on-a-chip, this con¬ 
cept has demonstrated precise single particle operations 
based on the selective motion of colloidal inclusions. 
Paramagnetic microspheres, which can be remotely con¬ 
trolled via non invasive magnetic fields, are currently em¬ 
ployed in biotechnological applications.The surfaces 
of these particles can be chemically functionalized allow¬ 
ing to bind selectively to defined targets.^ For an ensem¬ 
ble of monodisperse particles, the formation of a thresh¬ 
old energy where only a fraction of particles will move 
in response to an applied magnetic field is difficult to 
achieve. In most cases, the field-induced interactions be¬ 
tween the particles favor collective motion rather than 
selected displacements. Magnetic patterned substrates 
have shown considerable potential to overcome the above 
limitation. These patterns can generate strong localized 
field gradients, allowing controlled particle trapping and 
transport along predefined magnetic tracks. Recent ex¬ 
periments include the use of arrays of permalloy ellip¬ 
tical islands®, cobalt microcylinders,® domain wall con¬ 
duits,^® magnetic wires,exchange bias layer systems,^® 
and magnetic micromoulds. 

An alternative method consists in using epitaxially grown 
ferrite garnet films (FGFs), where magnetic domains 
with the width of few microns, i.e. on the particle 
scale, self-assemble into patterns of stripes or bubbles. 
Originally developed for magnetic memory applications^® 
and magneto-optical imaging,^® the FGFs are ideal to 
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manipulate paramagnetic colloids,or superconducting 
vortices.^® The highly localized driving force originates 
from the stray field gradient generated by the Bloch 
walls (BWs) in the FGF. When properly synthesized, 
the displacement of the BWs caused by an external field 
is smooth and reversible, with absence of wall pinning, 
thus creating a precise and controllable driving force for 
the particle motion. Here, we use a FGF to manipulate 
and transport an ensemble of paramagnetic particles or a 
fraction of it, demonstrating a technique to dynamically 
organize a colloidal system into trapped immobile parti¬ 
cles and particles which become mobile above a threshold 
field. 

A bismuth-substituted ferrite garnet of composition 
Y2.5Bio.5Fe5_qGagOi2 {q = 0.5— 1) was prepared by liq¬ 
uid phase epitaxial growth on a (lll)-oriented gadolin¬ 
ium gallium garnet (GGG) substrate. Oxide powders of 
the constituent elements, as well as PbO and B2O3, were 
initially melted at 1050 “C in a platinum crucible while 
the GGG wafer was located horizontally just above the 
melt surface. After lowering the temperature to 700“G, 
growth of the FGF was started by letting the substrate 
touch the melt. Keeping it there for 8 minutes pro¬ 
duced a FGF of 5 micron thickness, more details can 
be found in a previous work.^® At ambient temperature 
the FGF has a spontaneous magnetization perpendicular 
to the plane of the film. To minimize the magnetic en¬ 
ergy, the FGF breaks up into domains characterized by 
a labyrinth stripe pattern, easily observed by polarized 
light microscopy due to the large Faraday effect in this 
material. High frequency magnetic fields were used to 
transform this pattern into a regular triangular lattice of 
” magnetic bubbles”. These are cylindrical ferromagnetic 
domains magnetized oppositely to the remaining contin¬ 
uous area of the FGF, and separated by BWs. The mag¬ 
netic bubbles have in zero field a diameter of 6.4 /xm and 
a lattice constant of I = 8.6 ^m. Fig. 1(a). 

To avoid particle adhesion to the film, the FGF was 
coated with a 1 /xm thick polymer layer composed of 
a positive photoresist (AZ-1512, Microchem).®® Param- 
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FIG. 1. (a) Sketch of the magnetic bubble lattice (lattice con¬ 
stant I = 8.6/im) with paramagnetic colloids and subjected 
to a magnetic field rotating in the (a;; z) plane. Schematic on 
the right shows four magnetic bubbles (B) with two intersti¬ 
tials regions (/, 7) illustrating the two possible particle path¬ 
ways (dashed lines) and corresponding crystal directions, (b) 
Sequence of images showing the transport via particle swap¬ 
ping. (c) A particle defect transported via particle swapping, 
VideoSl. Superimposed are the trajectories of the particles, 
in green (blue) are trajectories along the BIB [BIB) path¬ 
way, VideoSl (Multimedia View), (d) Filling of a lattice va¬ 
cancy in the colloidal crystal, VideoS2 (Multimedia View). In 
these experiments, the fields have amplitudes = 0.7 kA/m, 
= l.QkAlm and frequency cu = 18.8s“^. 


agnetic microspheres with 2.8fj,m diameter (Dynabeads 
M-270) were diluted in highly deionized water (milli-Q, 
Millipore), and deposited above the FGF. After ^ 5 min 
of sedimentation, the particles became two-dimensionally 
confined above the FGF due to balance between repul¬ 
sive electrostatic interaction with the polymer layer and 
magnetic attraction. 

The magnetic field applied in the (x; z) plane was gener¬ 
ated by using two custom-made Helmholtz coils assem¬ 
bled above the stage of an upright microscope (Eclipse 
Ni, Nikon) equipped with a lOOx 1.3 NA objective and 
a 0.45X TV lens. The coils were arranged perpendicular 
to each other and connected to two independent power 
amplifiers (KEPCO BOP) driven by an arbitrary wave¬ 
form generator (TTi-TGA1244). 

In the absence of an external field, the BWs of the mag¬ 
netic bubbles attract the paramagnetic colloids and with¬ 
out the polymer coating, the particles sediment above 
the BWs. However, due to the polymer film, the parti- 
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FIG. 2. (a) Particle flux J versus normalized density p for dif¬ 
ferent frequencies, (b) J versus amplitude of the perpendicu¬ 
lar field Hz {Hx = 0.7 kAjm). Dashed red hues separate the 
regimes where all particles are immobile {Hz < 0.7 kAjm)^ 
only propagation of excess particles occurs (0.7 kA/m < Hz < 
1.2kA/m), and all particles are mobilized {Hz > 1.2kA/m). 
(c) J versus amplitude of the in-plane field Hx for three differ¬ 
ent values of Hz (d) Dependence of J on the driving frequency 
Lo for p = 1.27. 


cles have an higher elevation from the surface of the FGF 
and the magnetic potential become smoother, featuring 
energy minima at the centers of the magnetic bubbles. 
Thus, once deposited above the FGF, the particles form 
a perfect triangular lattice for a normalized areal density 
p = pa = 1, Fig.l. Here p = N/A is the particle number 
density, N is the number of particles located within the 
observation area A = 140 x 105 pm? and a = 64 pm? 
is the area of the Wigner Seitz (WS) unit cell around 
one bubble. For p> the excess particles redistributed 
within the magnetic domains, and each bubble became 
populated by colloidal doublets, triplets or larger clus¬ 
ters. 

We induce particle motion by applying an external 
magnetic field rotating in the {x]z) plane, H = 
(i7a; sin (wt), 0,77^ cos (w/:)), with angular frequency lo = 
2tlv and amplitudes {Hx;Hz), Fig.1(a). In most of the 
experiments, we keep fixed the amplitude of the in-plane 
component = 0.7 kA/m, and change the ellipticity of 
the applied field by varying the amplitude of the perpen¬ 
dicular component, Hz- 

For amplitudes 0.8 kA/m < < 1.2 kA/m, and load¬ 

ing p > 1, the particle transport takes place via a 
swapping mechanism in doubly occupied domains, where 
adjacent particles synchronously exchange their posi¬ 
tions. The bubble lattice thereby preserves the over¬ 
all occupancy of one particle per domain. Increasing 
the particle density, the swapping motion occurs in the 
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form of creation/destruction of doublets, triplets or even 
tetramers. However, for a wide range of particle densi¬ 
ties, we find that colloidal defects propagate mainly via 
doublets swapping motion, and the latter is illustrated in 
Fig. 1(b). 

To explain the mechanism leading to the defect propa¬ 
gation, let us consider the arrangement of four bubbles 
with their interstitial regions, as shown in right part of 
Fig. 1(a). Energy calculations^^ show that when the ro¬ 
tating field becomes anti-parallel with respect to the bub¬ 
ble magnetization, it generates in the interstitial regions 
two energy wells with triangular shape and opposite ori¬ 
entations, I and I, with corners pointing towards the 11 
and —1 — 1 directions, respectively. In this situation, 
there are two equivalent pathways along which an excess 
particle can propagate towards the —11 direction, either 
along the BIB pathway (dashed blue line), or along the 
BIB pathway (dashed green line). Both pathways are 
energetically equivalent, and the particle’s choice is dic¬ 
tated by the initial orientation of the doublet in the bub¬ 
ble domain. A doublet initially oriented along the —21 
direction will send a particle along the BIB pathway 
(Fig.1(b), first three images). Afterwards, this particle 
will form another doublet oriented along the —12 direc¬ 
tion, which will send a particle along the complementary 
BIB pathway (Fig.1(b), last three images). When the 
moving particle encounters a vacancy in the colloidal lat¬ 
tice, it will fill it and the defect propagation will end 
there, as shown in Fig. 1(c). 

We characterize the system conduction along the driv¬ 
ing direction by measuring the particle flux as J = p(u), 
where (v) is the average speed as determined from par¬ 
ticle tracking. Fig.2(a) shows J versus the dimensionless 
density p, for three different frequencies. Below the load¬ 
ing /5 = 1, i.e., having less than one particle per unit cell, 
there are no excess particles, and thus J = 0. For p > 1, 
we measure a net colloidal current which grows linearly 
with the loading up to p ~ 1.6. In this regime, only 
the excess particles contribute to the current, while one 
particle per magnetic bubble does not reach the mobility 
threshold. The speed of the particles inside the WS unit 
cell is phase-locked with the driving field, and given by 
V = Iuj/2tt where I is the lattice constant. Thus, increas¬ 
ing the driving frequency increases the average speed and, 
in turn, the slope of the curve in Fig.2(a). The direction 
of motion of the excess particles is dictated by the chi¬ 
rality of the rotating field, thus changing the polarity of 
one of the fields {Hx or H^) allows to invert the entire 
flow of particles across the film. 

Increasing the loading further (p > 2), the current 
reaches a maximum, and jamming between closely mov¬ 
ing colloids forbids further particle transport. Shown 
in Fig.2(b) is the effect that the perpendicular compo¬ 
nent of the applied field has on the particle current. 
This component acts directly on the size of the mag¬ 
netic domains, since the diameter of the magnetic bub¬ 
bles increases (decreases) when is parallel (antipar¬ 
allel) to the bubble magnetization. The graph, obtained 



FIG. 3. (a)-(f) Sequence of images showing the reversible 

assembly of a colloidal lattice by switching the applied field 
between the two threshold helds, and H^- First, the mag¬ 
netic particles were driven toward the 1 — 1 direction (green 
arrow in (b) by an applied field rotating in the (x; z) plane 
with components = Q.7 A/m and = lAA/m > _H|, 
and frequency a; = 18.8s“^. In (d) the direction of motion 
is inverted and the perpendicular component of the field de¬ 
creased to Hz = G.QkA/m < H!/- In (e) the held rotates 
in the (j/; z) plane, and the excess particles is transported to¬ 
wards the — 1 — 1 direction, Filing the whole lattice of bubbles, 
VideoSS (Multimedia View). 


for a fixed density of p = 1.27, shows that the parti¬ 
cle flux grows in discrete steps as the amplitude of Hz 
increases. Below Hz = H^ = 0.7 kA/m, no current is 
observed. Increasing Hz the excess particles start to be 
mobilized, and the flux raises till reaching the constant 
value J = 6.1 for Hz > 0.7kA/m. Increasing 

Hz further, reveals a second field value, = 1.2 kA/m, 
where the periodic displacements of the BWs are able 
to drive all particles synchronously across the film. The 
field values Hf and are therefore the mobility edges 
where different sub-ensembles of particles can be set into 
motion. 

The effect of the in-plane component of the applied field 
Hx on the flux J is shown in Fig. 2(c). While Hz con¬ 
trols the size of the magnetic domains, the effect of Hx is 
to break the spatial symmetry of the potential, inducing 
a net particle current towards a defined direction. When 
Hz < HI, no current is observed for any value of Hx. For 
Hz = 1.0 kA/m the flux is composed of excess particles, 
while for Hz = 1.3kA/m all particles are mobilized and. 
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in both cases, when > 0.7 kA/m it becomes indepen¬ 
dent on further increasing of Hx- 

Fig.2(d) shows the dependence of J on the driving 
frequency for p = 1.27. The current displays a 
full participation of the excess particles up to w ~ 
44 where the flux reaches its maximum value, J = 
13.8 corresponding to an average defect speed 

of ri = 60 p,ms~^. Beyond this frequency, the participa¬ 
tion of excess particles becomes partial, and the current 
decreases monotonously, reaching zero near w = 110 s“^. 
In this second regime, the high frequency motion of the 
particles is found to be intermittent, with the excess par¬ 
ticles randomly switching between being immobile for 
some period, for afterwards becoming fully mobile again, 
reducing the efficiency of our magnetic device. 
Controlling the motion of the excess particles makes it 
possible to easily create or destroy a colloidal lattice 
by switching the applied field between the two mobil¬ 
ity thresholds. This concept is demonstrated in Fig.3 
(Videos). Starting from a triangular arrangement with 
one particle per WS cell (a) {t = Os), we apply an ellip- 
tically polarized magnetic field with components = 
0.7 kA/m, Hz = 1.3 kA/m > i7|, and angular frequency 
w = 18.8s“^, which mobilize the whole lattice at a speed 
of 25.8p,ms~^ towards the 1 — 1 direction, Figs.3(b,c), 
leaving the bubble substrate almost unfilled. Since the 
moving particles are phase-locked with the driving field 
for the used frequency, the translating lattice is stable 
and preserves the initial triangular order during motion. 
After t ~ 10s, we change the polarity of the in-plane field 
{Hx —> —Hx), in order to invert the particle flux towards 
the —11 direction. We also decrease the held amplitude 
to Hz = O.OkA/m < i7|, inducing a current composed 
only of excess particles which start filling again the bub¬ 
ble lattice in the bottom part of the film, Figs.3(d,e). 
Since, during these operations the top part of the film is 
left unfilled, we change the orientation of the applied field 
after t = 26.7s. The field now rotates in the {y, z) plane, 
and transports the excess particles toward the —1 — 1 
direction, reforming the colloidal crystal after 31s. The 
colloidal assembly demonstrated here can be further opti¬ 
mized by either controlling the particle density, and thus 
the amount of excess particles propagating through the 
lattice, or by increasing the driving frequency, up to a 
maximum speed of u = 60 p,ms~^, which will further 
speed up the re-writing process. 

In conclusion, we demonstrate a technique to remotely 
generate and control the motion of defects in two dimen¬ 
sional lattices, while keeping track of the position of the 
individual particles, which could be used as model system 
to study the dynamics of impurities in crystalline materi¬ 
als. The energy scales involved in our system (^ IbOksT, 
with T ~ 293K) are much beyond the effect of ther¬ 


mal fluctuations, which could interfere with the colloidal 
transport process. This makes our magnetic device fully 
controllable, ensuring an extremely precise tuning of the 
particle speed and dynamics in real time and space. Al¬ 
though our experiments focus on using FGF films as func¬ 
tional platform, the technique reported here should be 
applicable, within the constraint of colloidal particle size 
and lattice wavelength, to other platforms where mag¬ 
netic patterns are created by ’’top-down” fabrication pro¬ 
cesses. 
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